Predicted Ultrafast Single Qubit Operations in Semiconductor Quantum Dots 
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Several recently proposed implementations of scalable quantum computation rely on the ability to 
manipulate the spin polarization of individual electrons in semiconductors. The most rapid single- 
spin-manipulation technique to date relies on the generation of an effective magnetic field via a 
spin-sensitive optical Stark effect. This approach has been used to split spin states in colloidal CdSe 
quantum dots and to manipulate ensembles of spins in ZnMnSe quantum wells with femtosecond 
optical pulses. Here we report that the process will produce a coherent rotation of spin in quantum 
dots containing a single electron. The calculated magnitude of the effective magnetic field depends 
on the dot bandgap and the strain. We predict that in InAs/InP dots, for reasonable experimental 
parameters, the magnitude of the rotation is sufficient and the intrinsic error is low enough for them 
to serve as elements of a quantum dot based quantum computer. 

PACS numbers: 03.67.Lx, 73.63.Kv 



Scalable proposals for quantum computation [0 require 
well-defined individual qubits that can be manipulated 
individually (single-qubit gates) and also can control- 
lably interact with each other (two-qubit gates) ||. In 
solid-state systems often the two-qubit gates appear eas- 
ier to realize, e.g. by the exchange interaction between 
electrons The manipulation of single qubits has been 
perceived as more difficult; this motivated the proposal 
of all-exchange-based quantum computation [Q, which re- 
quires only two-qubit (exchange) gates. All-exchange- 
based computations require a large number of gate oper- 
ations, although recent work indicates substantial reduc- 
tions are possible Q." 

Here we predict that a spin- AC Stark effect (1) pro- 
duces coherent rotations of electron spin in quantum 
dots, (2) the rotation angle can exceed tt for reasonable 
experimental parameters, and (3) the error rates are tol- 
erable for quantum computation. The spin-splittings of 
states in CdSe quantum dots||, j?j and coherent manipu- 
lation of ensembles of spins in ZnMnSe quantum wells j^] 
has been demonstrated experimentally. Theoretical con- 
siderations of the AC Stark effect || [II], |ll| have focused 
to date only on non-spin-selective shifts of energy lev- 
els. Here we consider a spin-selective AC Stark effect 
whereby a quantum dot is illuminated with a single in- 
tense pulse of circularly-polarized nonresonant light (Fig. 
la). Such a pulse shifts the energies of dot states, and 
due to differing transition matrix elements the two spin 
states are shifted differently. Hence this pulse produces a 
splitting of the two lowest energy conduction states. The 
splitting of these two states (spin-up and spin-down) can 
be viewed as an optically-induced pseudo-magnetic field 
( B e ff ) acting on the state pair. Direct application of 
this approach to quantum computing is clear, for one 
well-known physical realization of a single qubit opera- 
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FIG. 1: a) SAQD and CQD quantum dot geometries. For 
SAQDs, the incident light is taken to be incident along the 
growth direction, b) Level diagram for a quantum dot with 
one electron in the lowest conduction band state, perturbed 
by incident light detuned below the bandgap. The electron 
spin precesses in response the the effective magnetic field B e f f 
corresponding to the spin splitting. 



tion is a magnetic field applied to a spin for a definite 
period of time. We find that the size of the effective field 
in certain dot systems, for reasonable experimental con- 
ditions, exceeds the magnitude required for 180 degree 
rotation of the spin polarization (7r-pulse). Furthermore 
we evaluate the error rates for spin manipulation via this 
process under reasonable conditions to be smaller than 
10 -6 , which is within the tolerance for quantum compu- 
tation. 

There are a wide variety of quantum dot systems from 
which to choose, with various properties affecting the size 
of the Stark splitting. Of particular interest are the level 
spacings and oscillator strengths, both of which are af- 
fected by material parameters and the size of the dots. 
We will examine Stark splittings in CdSe and InAs col- 
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Parameter 


InAs 


GaAs 


InP 


CdSe 


E g (eV) 


0.418 


1.519 


1.424 


1.84 


A so (eV) 


0.38 


0.341 


0.11 


0.42 


VBO (eV) 


0.085 


0.0 


-0.045 


_ 


* 

m e 


_ 


_ 


_ 


0.119 


71 


19.67 


6.98 


4.95 


2.52 


72 


8.37 


2.25 


1.65 


0.65 


73 


9.29 


2.88 


2.35 


0.95 


E P (eV) 


22.2 


22.7 


20.4 


17.4 


a c (eV) 


-6.67 


-9.55 


-7.33 


_ 


a v (eV) 


1.67 


0.95 


0.73 


_ 


6 (eV) 


-1.8 


-2.0 


-2.0 


- 


a (eV) 


-3.6 


-5.4 


-5.0 




az a t t (nm) 


0.6058 


0.5653 


0.5869 




Cxxxx (GPa) 


832.9 


1211 


1022 




Cxxyy (GPa) 


452.6 


548 


576 




Cxyxy (GPa) 


395.9 


604 


460 





TABLE I: Material Parameters used in the calculations[[L2|. 
The valence band offset (VBO) is the unstrained valence en- 
ergy, determining the band alignment in SAQDs. For CQDs 
infinite barriers were used for the vacuum. 



the larger bandgap and lack of strain, the CdSe cal- 
culations were performed with a single-band model for 
the conduction band, and a four-band model for the va- 
lence band. The CdSe calculations were done for the 
zincblende form. The material parameters are given in 
Table I. 

Calculation of the energy shifts under illumination was 
done by constructing a restricted set of Fock states for the 
combined electron-photon system. Eight valence states 
and the lowest conduction doublet were used for the 
quantum dot states. The states in the basis for the 
electron-photon system were |1,0;1, 1,1, 1,1, 1,1,1; 7V 7 ), 
10,1,1,1,1,1,1,1,1,1;^), |1,1;1,1,...0,...,(JV-1) 7 > 
for a total of 10 states. For the first two Fock states all 
valence states and one conduction state are occupied, and 
there are 7V 7 photons in the photon field. For the other 
eight Fock states, both conduction states are occupied, 
only one valence state is unoccupied, and there is one 
less photon than in the first two states. Inclusion of ad- 
ditional states did not alter the results. In this restricted 
basis the Hamiltonian for the coupled electron-photon 
system is given by 



loidal quantum dots (CQD), as well as in InAs/GaAs and 
InAs/InP self assembled quantum dots (SAQD). This set, 
while not exhaustive, spans a broad range of parameters 
by including unstrained spherical and strained lens-like 
dots with widely varying bandgaps. 

The AC Stark effect with unpolarized light is a nonlin- 
ear effect whereby light with photon energy tuned near 
to the absorption transition between two states induces a 
level repulsion between the two states^, [H]. In pertur- 
bation theory [|Io[ this level repulsion depends on If/S, 
where / is the light intensity, / is the oscillator strength 
of the transition, and 5 is the detuning of the photon en- 
ergy from the transition energy (see Fig. lb). When 
circularly polarized light illuminates a transition from 
the first valence state pair to the first conduction state 
pair in a crystal with spin-orbit interaction the oscilla- 
tor strengths of the transitions differ considerably; this is 
the source of the effective spin splitting of the conduction 
state pair. 

Optical Stark shifts were calculated non-perturbatively 
using a restricted basis of quantum dot wave functions 
calculated with k • p theory in the envelope approxima- 
tion. The wave functions were calculated numerically on 
a real-space cubic grid using the Lanczos algorithm [[L3|. 
This method has been used previously to calculate a va- 
riety of electronic properties of strained QDs. The grid 
spacing was chosen to match the lattice constant of the 
barrier material in the SAQDs, and the dot material it- 
self for the CQDs. The number of grid sites varied from 
10 4 -2x 10 5 depending on the size of the QD. For SAQDs 
the strain arising from lattice mismatch was calculated 
by doing a conjugate gradient minimization of the con- 
tinuum elasticity strain energy. The strain was then used 
as input to an eight-band strain-dependent Hamiltonian, 
which was diagonalized as described above. Because of 
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where P is the momentum operator, e is the polarization 
vector of the incident light with electric field Eq, and 
^cj are the valence and conduction band states respec- 
tively, and d^ are the dipole matrix elements. Coherent 
states for the light field can then be naturally constructed 
from the mixed electron-photon states obtained above. 

In Fig 2a and 2b are shown the energies of the con- 
duction states for two representative dot systems: an 
InAs/GaAs SAQD and a CdSe CQD. In both cases,the 
splittings increase as the detuning approaches zero, and 
change sign when the detuning changes sign. The first 
notable feature is the difference in spin selectivity, which 
can be traced to the differing strain and shape. For CdSe 
both spin states are shifted, but by different amounts, 
while for InAs/GaAs only one spin state is shifted. Be- 
cause the CdSe dot has approximate spherical symmetry 
and is unstrained, the highest valence states are closely 
spaced ( < 1 meV), with comparable amounts of heavy 
hole (HH) and light hole (LH) character. As a result, 
both the HH and LH transitions contribute to the Stark 
shift, and both spin directions for the conduction state 
are shifted. The difference between the shifts reflects the 
difference in oscillator strengths for HH and LH transi- 
tions, whose ratio of 2:1 is approximately the ratio found 
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FIG. 2: a) Energy levels for an InAs/GaAs dot with 1.245 eV 
bandgap. The incident power density is 10 9 W/cm 2 and the 
conduction band electron has spin parallel to the direction of 
the incident light, b) CdSe dot with 2.030 eV bandgap. 



in bulk band-to-band transitions, 3:1. In contrast, for 
InAs/GaAs dots the highest valence state is almost en- 
tirely HH, and separated from the next valence state by 
several tens of meV (depending on size). Hence, the Stark 
shift is dominated by the (doubly degenerate) highest va- 
lence state which gives a Stark shift for only one spin 
direction of the conduction electron. 

We now examine the effect of dot size on the spin split- 
ting. Shown in Fig. 3 are the splittings for detunings 
from 30 meV - 70 meV for the four systems under consid- 
eration. Since we are ultimately interested in manipulat- 
ing spins through the effective magnetic field, it is useful 
to consider the precession angle associated with a light 
pulse of duration St , given by S = AESt/h where AE is 
the Stark spin splitting between up and down conduction 
states. In the results that follow, we give Stark splittings 
in both meV, and the corresponding S for a 200 fs pulse 
with power density 10 9 W/cm 2 , which we refer to as a 
reference pulse. The ~ 2meV spin splitting seen in CdSe 
CQD's agrees with that measured experimentally [||. 

The trend in Fig. 3 is surprising; for a given material 
system larger dots have larger splittings, in spite of the 
fact that oscillator strengths decrease with increasing dot 
size. The reason is that the the interaction term in the 
Hamiltonian is proportional to 1/uj. For a fixed detuning, 
the increase in 1/uo for larger dots dominates the decrease 
in the dipole matrix element. 

The importance of the bandgap suggests selecting the 
dot material to minimize the gap. One possibility is to 




FIG. 3: Stark splittings as a function bandgap and detuning 
for InAs, InAs/InP, InAs/GaAs, and CdSe quantum dots. 
The angle S = AESt/h is for an incident 200 fs pulse with 
10 9 W/cm 2 . 
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FIG. 4: Sensitivity of S to variations in the bandgap as a 
function of bandgap and detuning, estimated using a finite 
difference approximation. Note that larger detunings than in 
Fig. 2 were used so there would be values of S at fixed values 
of hoo for every bandgap value. S is computed assuming a 
reference pulse 



use a small bandgap material, such as InAs, but to avoid 
the bandgap-increasing effect of strain by using an InAs 
CQD. Confinement increases the bandgap over that of 
bulk InAs, but for typical CQD sizes E g < 1 eV, which 
is substantially less than E g w 2 eV for CdSe CQDs. 

Fig. 3 confirms the expectation that InAs CQDs have 
larger splittings than CdSe CQDs over most of the range 
of sizes, but the increase in splitting is at most a factor of 
1.5. While InAs CQDs have larger splittings than CdSe, 
they are still a factor of 5 smaller than for InAs/GaAs. 
In spite of a smaller bandgap, the lack of strain in InAs 
CQDs decreases the HH/LH splitting so much that the 
Stark splitting is substantially decreased as well. 

An alternative to the complete elimination of strain 
is to choose a substrate material that has a lattice con- 
stant closer to InAs. InAs/InP is an excellent candi- 
date, with a lattice mismatch of 4% (as opposed to 7% 
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for InAs/GaAs), and a edge-to-edge bandgap of approx- 
imately 0.5 eVj[L4|, [l^l . The highest valence state is still 
predominantly HH, and the first excited valence state 
ranges from 10 meV - 30 meV below, depending on size. 
In addition, InAs/InP has a bandgap at the technolog- 
ically important 1.55 /im. Fig. 3 shows that the split- 
tings in InAs/InP are substantially larger than those for 
InAs/GaAs. Some of the improvement is simply due to 
the smaller bandgap and the factor of 1/oj in the inter- 
action as discussed earlier. However, the discontinuity in 
the curves indicates the improvement due to the increase 
in strain and the resulting increased HH/LH splitting. 
At our standard illumination of 10 9 W/cm 2 with a 200 
fs pulse width, a 7r-rotation is possible at even large de- 
tunings, 8 = — 70meV. 

We now consider errors in S due to the (unavoidable) 
variations in the bandgap of the dots. For measurements 
on an ensemble of dots, inhomogeneities in the dot size 
will yield a different rotation angle for different dots. 
More relevant for a quantum computer, in which indi- 
vidual dots will be selected, the finite line width of the 
states will cause some uncertainty in S . To address this 
question, we have calculated numerical estimates of 

1 dAE(E g ,hm) = 1 d0 s (E g ,hm) 
AE dE g S dE g U 

as shown in Fig. 4. For the materials and sizes 
considered this quantity varies from —0.008 meV -1 to 
—0.014 meV -1 , depending on detuning. The results of 
Fig. 4 may be used to estimate the uncertainty in the 
rotation angle, A0 S . For example, an ensemble measure- 
ment of InAs/InP dots with an inhomogeneous linewidth 
of 50 meV would give A0 S ^0.5 for a 7r-pulse, which 
should be sufficient to permit observation of a spin echo. 
For a single dot with a linewidth of 0.1 meV, A6 S ^0.001, 
corresponding to a bit error rate of 10 -6 . This is suffi- 
ciently small for error correction algorithms to apply (l6| 
(desirable error thresholds have been estimated |l(| at 

10~ 5 - io- 6 ). 

Errors in S may also arise from shot noise in the laser 



pulse. The electric field E Q appearing in Eq. 2 has some 
uncertainty due to variations in the number of photons in 
the pulse. Assuming the laser is focused to 1 fim 2 with a 
photon energy of 1 eV, a reference pulse contains approx- 
imately 10 T photons. The uncertainty in the incident in- 
tensity is AI « I/y/N 7 , giving A6 S « 2tt/VN 7 w 10~ 3 . 
The corresponding bit error rate is approximately 10 -6 , 
which is still acceptably small for error correction [|l6|. 
This estimate gives a lower limit on the error rate since 
it neglects additional sources of laser noise that would 
increase the error rate. 

We conclude that use of the spin- AC Stark effect is a 
viable approach to single qubit manipulation in quantum 
dots. The magnitude of the Stark splitting in SAQDs is 
5-10 times larger than in CQDs due to the strain-induced 
HH/LH splitting. The Stark splitting also increases with 
decreasing bandgap, though the effect is smaller that of 
strain. Because strain increases both the HH/LH split- 
ting and the bandgap, finding the optimal system in- 
volves a seeking a system with sufficient strain to in- 
duce HH/LH splitting, but not so much as to make the 
bandgap too large. We propose as a candidate system 
InA/InP dots, which are strained, but have a relatively 
small bandgap of 0.7-0.9 eV. This energy range includes 
1.55 /im, which may be important for interconnection 
within a quantum computer, or for quantum commu- 
nication applications. We find that for such dots, tt- 
pulses may be obtained for experimentally realistic pulses 
(10 9 W/cm 2 , 200 fs) and detunings (-70 meV). For typ- 
ical inhomogeneous broadening (50 meV) the variation 
in rotation angle is A0 S « 0.5. For a single dot with an 
intrinsic linewidth of 0.1 meV we estimate the bit error 
rate to be on the order of 10 -6 , and we estimate the bit 
error rate due to laser shot noise is also on the order of 
10~ 6 . These error rates are sufficiently low for quantum 
error correction. 
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